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ABSTRACT 

We review the physical processes occurring in the magnetosphere of accreting X-ray pulsars, 
^vj ' with emphasis on those processes that give rise to observable effects in their high (E>10 keV) 

energy spectra. In the second part we compare the empirical spectral laws used to fit the observed 
spectra with theoretical models, at the light of the BeppoSAX results on the broad-band charac- 
terization of the X— ray pulsar continuum, and the discovery of new (multiple) cyclotron resonance 
features. 
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1. INTRODUCTION 

43 
CLi An X-ray pulsar is, by definition, a celestial source showing pulsed emission when 

observed in X-rays. The very first pulsating X-ray source, Centaurus X-3, was 

discovered in 1971 by the first scientific X-ray satellite Uhuru ||. Its 4.8 s pulse 

c/3 . period implied a small emitting region, and because the object responsible for the 

pulsation is not destroyed by the centrifugal force it is necessary that at its surface 

the gravitational force is greater than the centrifugal one. This implies Sl p ^ y/Gp, 

where Q p is the pulse frequency, G the gravitational constant, and p the object 

mean density. The observed value of £l p implies p <; 10 7 g/cm 3 and therefore the 

compact nature of the object responsible of the pulsed emission was established. 

The binary nature of Cen X-3 was soon after recognised by the observation of 

Doppler modulation in the observed pulse period |2q ]. The 2.1 day modulation was 

coincident with the periodic disappearing of the source X-ray flux, interpreted as 

eclipse of the compact object by the companion. Finally the optical counterpart was 

discovered as an early-type O star |Q. With all these elements it was possible to 

determine the mass of the compact object, which resulted to be 1.4 M : a neutron 

star (NS). 
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2. PHYSICAL PROCESSES IN X-RAY PULSARS 

The physical scenario able to explain the production of pulsed X-ray emission was 
elaborated by Shklovskii [E7[ before the discovery of Cen X-3. X-rays are produced 
in the conversion of the kinetic energy of the accreted matter (coming from the 
intense stellar wind of an early-type star — wind-fed binaries, or coming from 
an accretion disc due to Roche-lobe overflow — disk-fed binaries) into radiation, 
because of the interaction with the strong magnetic field of the NS, of the order of 
10 11 -10 13 gaussR. The dipolar magnetic field of the NS drives the accreted matter 
onto the magnetic polar caps, and if the magnetic field axis is not aligned with the 
spin axis, the NS acts as a "lighthouse", giving rise to pulsed emission when the 
beam (or the beams, according to the geometry) crosses our line of sight. 

For a detailed description of the spectral properties of accreting X-ray pulsars 
(AXPs) it is therefore necessary to describe the interactions of the X-rays produced 
at the NS surface with the highly magnetized plasma forming the magnetosphere. 
This is a formidable task because we cannot use a linearized theory for the radia- 
tive transfer equations but we have to deal with the fully magnetohydrodynamical 
system. This is due to the fact that the coupling constants among the interactions 
are so large that a series expansion is impossible. 

This is the reason why there is not a parametrized description of AXP spectra in 
terms of physical quantities, but only empirical laws to fit the observed spectra. An 
alternative method is the numerical solution of the radiative transport equations as- 
signing particular values to the physical parameters, comparing the obtained spectra 
with the observed ones, and varying the parameters until a match is reached. 

2.1. Cyclotron resonant features 

At some distance from the NS, that we will call magnetospheric radius r m , the 
motion of the accreted matter will be dominated by its intense magnetic field. We 
define magnetosphere the region around the NS delimited by r m . The electrons 
present in the magnetosphere will have an helicoidal motion along the magnetic 
field lines, with gyromagnetic (Larmor) frequency given by 

(i) 



where 7 is the Lorentz factor. For the magnetic field strength B expected in the NS 
magnetosphere, the motion of the electron in the direction perpendicular to B is 
quantized in the so-called Landau levels (see e.g. p3[ ). In the nonrelativistic case, 
the energy associated to each level is given by 

fiuj n — n hto c (2) 



1 Obtained from conservation of magnetic flux during the process of collapse from a "normal" 
star (B ~ 10-100 gauss, R ~ 10 6 Km) to a NS (R ~ 10 Km) 



where u> c is the Larmor gyrofrequency given in Eq. [0. As an aside, from Eq. 
we have that E n = 11.6 • -B12 keV, where -B12 is the magnetic field strength in 
units of 10 12 gauss. Therefore we expect to observe cyclotron features in the hard 
(E > 10 keV) energy range. As we have seen, in the non-relativistic case the energy 
levels are harmonically spaced. When relativistic corrections are taken into account 
a slight anharmonicity is introduced in the Landau levels. Indeed, we have 

2 V'mc' 2 + 2nhuj c sin 9—1 

nu) n = mc „■ (3) 

sin 9 

where 9 is the angle between the line of sight and B. 

Another consequence of the existence of the Landau levels is that an electro- 
magnetic wave propagating in such a plasma will have well defined polarization 
normal modes, i.e. the medium will be birifringent m. It is not our intention to 
enter into the details of the propagation of waves in the magnetospheric plasma; 
we will develop a semi-quantitative approach by highlighting the plasma properties 
that have observable consequences in the AXP spectra. 

If we introduce the complex refraction index N, with its real part the geometric 
refraction index and with its imaginary part the absorption coefficient, then the 
dispersion relation in the non-relativistic case can be written as a bi-quadratic 
equation in N. The solution for N will have the form JL| 

N? ex — !— NS oc — !— . (4) 

CiJ — CiJ c U> + UJ C 

The wave with Ni presents resonance and is right-handed circularly polarized 
(that is in the same sense as the electron gyration). This wave is called extraordi- 
nary, in opposition to the ordinary wave described by N2, which is left-handed 
circularly polarized. By introducing the complex refractive index is straightforward 
to obtain the cyclotron absorption cross section. By means of the optical theorem 
we obtain 

a c = Air 2 af — \ei\ 2 S{uj -lu c ) (5) 

where a/ = e 2 /he is the fine structure constant, and el is the polarization versor 
of the extraordinary wave. 

Up to now, we worked neglecting both relativistic corrections and thermal mo- 
tions (cold plasma approximation). The release of the latter condition allows an 
electron to absorb waves not only of frequency u> = w c , but in the interval uj c ± Alud, 
where the Doppler width is given by 



2kT 

Au>d=Uc\ 9 I cos 6»| (6) 

V mc z 

where T is the electron temperature (we assumed a Maxwell-Boltzmann distribution 
for the electrons). 



Once the electron absorbs a photon it (almost) immediately de-excitates on a 
time scale t r ~ 2.6 x 10 -16 B^ sec jig ]. This has important consequences for the 
scattering cross sections. Indeed, while a scattering process involves two photons 
(one going in, one going out), absorption (or emission) processes involve only one 
photon. Therefore one expect that the two cross section are different. This is not 
true just because an absorbed photon is immediately re-emitted, and therefore the 
absorption-emission process is equivalent to a scattering. It is possible to show ]13| 
that the cyclotron scattering cross section cr res has the same form as the cyclotron 
absorption cross section a c (Eq. |5|) with the prescription 

6{uJ " " c) -* (u-^y + rt/4 (7) 

where r r = "f r tu c , and ~f r = (2/3) (e 2 /mc 3 )u) is the radiative damping. 

Therefore photons with frequency close to lo c will be scattered out of the line 
of sight, creating a drop in their number. Cyclotron "lines" observed in the spectra 
of AXPs are therefore not due to absorption processes, but are due to scattering of 
photons resonant with the magnctosphcric electrons (as it occurs for the Fraunhofcr 
lines in the Solar spectrum). This is why we will not use the term cyclotron lines 
but the more appropriate "cyclotron resonant features" (CRFs). 

When relativistic effects are taken into account, it is possible to show that also 
ordinary waves show resonance, and are scattered out of the line of sight. Another 
important effect on the radiative properties of the plasma is due to a pure quantum 
effect: the so-called vacuum polarization. The magnetic field at which a quantum 
mechanical treatment of the plasma is necessary can be defined when the classical 
cyclotron energy hu> c becomes equal to the electron rest mass mc 2 . That is 

Tn c 
B cr = — — = 4.414 x 10 13 gauss . (8) 

We call B cr the "critical" magnetic field strength. For B not far from B cr virtual 
electron-positron pairs can be created. These virtual photons dominate the polar- 
ization properties of the plasma for frequencies in the range u v \ Si u> Ss u> V 2, where 

id ra 



ui vX ~ 3keV ; lu v2 ~u} c (9) 

±*/U.li5 cr 

and therefore affect the scattering cross sections in-ii is the electron density in units 
of 10 22 cm -3 ). In Fig. [Owe show the effects of the inclusion of vacuum polarization 
on the opacity (cross section times density) as obtained by J32| . 

2.2. Continuum emission 

The main physical process responsible for the continuum emission in AXPs is Comp- 
ton scattering. We will not enter into the details of the problem of repeated scat- 
terings in a finite, thermal medium (see e.g. pq|). Let us only summarize that an 
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FIGURE 1. Angle-averaged scattering opacities (cross section times density) for 
the extraordinary (p\) and ordinary (172) modes, computed for 71,22 = 1 and %ui c = 
50 keV. Left: without vacuum polarization; Right: with vacuum polarization. From 
Ventura et al. (1979). 



input photon of energy E^ will emerge from a cloud of non-relativistic electrons (at 
a temperature T) with an average energy Et ~ E t e v (this is valid in the regime 
Ef <C 4fcT). The comptonization parameter y therefore gives a measure of the 
photon energy variation in traversing the plasma, and is given by 



4kT 



mc^ 



max (t, t ) 



/AkTY 



- I max (t, t 
\mc J 



Nonrelativistic 



Relativistic 



(10) 



where max (r, t 2 ) is nothing else but the average number of scattering suffered by 
the photons (t is the optical depth of the medium). Note that if J5j < 4fcT then 
photons can increase their energy at the expense of the electrons: this is inverse 
Compton scattering. 

The detailed description of the spectrum of the emergent photons requires the 
solution of the Kompaneets equation, but it is possible to obtain qualitative infor- 
mation for special cases: 

• y <C 1 In this case only coherent scattering is important, and the emergent 
spectrum will be a blackbody spectrum or a "modified" blackbody spectrum 
according whether the photon frequency is lower or greater than the frequency 
at which scattering and absorption coefficients are equal E3]. 
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FIGURE 2. Spectrum due Compton scattering in a thermal, nonrelativistic medium. 
Left: At low frequency the spectrum is blackbody, while at higher frequencies de- 
velops a Wien hump due to saturated inverse Compton. Right: Spectrum produced 
by unsaturated inverse Compton scattering (adapted from |23|j). 



• y ~S> 1 Inverse Compton scattering can be important. If we define a frequency 
u> co such that y(co co ) — 1, then for u) ^ uj co the inverse Compton scattering 
is saturated and the emergent spectrum will show a Wien hump, due to low- 
energy photons up-scattered up to hui ~ 3kT [£3| . In the case in which there 
is not saturation a detailed analysis of the Kompaneets equation shows that 
the spectrum will have the form of a power law modified by a high energy 
cutoff [23, 29]. These two regimes are qualitatively depicted in Fig. @. 



3. SPECTRAL X-RAY OBSERVATIONS OF AXPS 

3.1. Before BeppoSAX 

The very first observation of a CRF in a spectrum of an X-ray pulsar was per- 
formed in 1978 when Triimper et al. pl| observed a ~35 keV CRF in the spectrum 
of Hercules X-l. A while later it was observed not only the fundamental but also 
the first harmonics in the spectrum of the transient X-ray pulsar 4U0115+63 ]33] ]. 
Observations of CRFs in other AXPs showed that they are a quite common phe- 
nomenon in this class of objects. But it was with the advent of the Japanese satellite 
Ginga that a systematic analysis of the spectra of AXPs was performed in search 
of CRFs. Mihara JlfJ analysed the spectra of 23 AXPs and found that 11 among 
them showed CRFs. 



3.1.1. Continuum characterization 

Because CRFs are broad features, the exact determination of the continuum is 
of paramount importance. From the analysis of the HEAO-1/A2 spectra of AXPs 
White et al. |34| found an empirical law that was able to fit their energy spectra 



E E < E cut 

VOm { E) = { ^^E-E^ E>Ecut (11) 

It is evident that this model tries to simulate the unsaturated inverse Comp- 
ton process shown in Fig. 0. But this model suffers the problem of a too abrupt 
break around the cutoff energy E cut , problem enhanced by following more sensitive 
instruments. Therefore Tanaka p(| introduced a "smoother" cutoff of the form 

FBCO(E) = Je=E c . n (12) 

1 + exp 



Ef 

that he called Fcrmi-Dirac cutoff because of it resemblance with the Fermi-Dirac 
distribution function. It is important to stress that the FDCO model does not have 
any physical meaning: it only gives a better description of the break in the AXP 
spectra. 

Makishima and Mihara |ll]] were the first to note that in the AXPs showing 
CRFs there was a correlation between the cutoff energy E cut and the CRF energy 
E c , namely E c ~ (1.2 — 2.5) ■ E cut . Therefore it seemed that the cutoff was in some 
way due to the presence of the CRF. The next step was performed by Mihara, who 
introduced the so-called NPEX (Negative Positive Exponential) model 

NPEX(£) = (AE- a + BE +P ) exp ( —^- J . (13) 

This model is quite successful in describing the AXP spectra observed by Ginga 
in the 3-30 keV. Its components have also a physical meaning, because it mimics the 
saturated inverse Compton spectrum shown in Fig. if (3 = 2. Furthermore, because 
the (non relativistic) energy variation of a photon during Compton scattering is |23j 

AE _AkT-E 

~E~ ~ mc 2 (1 ' 

then when E — E c the medium is optically thick and therefore E c ~ AkT . 

3.1.2. CRF characterization 

Mihara, besides the introduction of the NPEX model to describe the AXP contin- 
uum, introduced a new form for the CRF 

which has the form of a Lorenzian of width W , and depth r. From a physical point 
of view, this is the form assumed by the cyclotron scattering cross section described 
by Eq. 



Source 


Obs Date 


E cyc (keV) 


FWHM (keV) 


References 


4U0115+63 (M) 


20 Mar 1999 


12.78 ±0.08 


3.58 ±0.33 


B 


4U1538 52 (M) 


29 Jul 1998 


21.5 ±0.4 


6.7 ±1.2 


|22] 


Cen X-3 (M?) 


27 Feb 1997 


28.5 ±0.5 


7.3 ±1.9 


\24\ 


XTE J1946+27 


09 Oct 1998 


33 ±4 


16 ±2 


|2£| 


OA01657 415 


04 Sep 1998 


36 ±2 
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<} 


4U1626 67 


06 Aug 1996 


38.0 ±0.9 


11.8 ±1.7 
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4U1907+09 (M) 


29 Sep 1997 


38.3 ±0.7 


9.7 ±2.3 


|1| 


Her XI 


27 Jul 1996 


42.1 ±0.3 


14.7 ±1.1 


ISL 


GX301 2 (M) 


24 Jan 1998 


49.5 ±1.0 


17.9 ±2.5 


|2C] 


Vela XI (M) 


14 Jul 1996 


54.8 ±0.9 


25.0 ±2.1 


|1£] 


GX1+4 


25 Mar 1997 










GS1843+00 


04 Apr 1997 
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X Persei 


09 Sep 1996 











M stands for multiple lines detected/suspected 

TABLE 1. BeppoSAX observations of X-ray binary pulsars 

3.2. BeppoSAX observations of AXPs 

With the advent of BeppoSAX the study of energy spectra of AXPs received new 
impulse because it was now possible to characterize with unprecedent detail the 
continuum on a broader energy range (0.1-200 keV), and the two high energy in- 
struments aboard BeppoSAX, namely HPGSPC (sensitive in 5-60 keV; |T^]) and 
PDS (15-200 keV; |jj), are the best suited for the detailed spectroscopy of CRFs. 
BeppoSAX observed all the persistent AXPs, plus a couple of transient ones (see 
Table fl]). As a first result, we found that the NPEX model, successfully used to 
fit the Ginga data, is not adequate to describe the broad AXP continuum ||. In 
particular we find that their continuum can be described in terms of (i) a black- 
body component with temperature of few hundreds eV; (ii) a power law of photon 
index ~1 up to ~10 keV; and a (hi) a high energy (>10 keV) cutoff that makes the 
spectrum rapidly drop above ~40-50 keV. 

Furthermore, the CRFs observed with BeppoSAX are better described in terms 
a Gaussian in absorption, defined as P8fl 



GAUABS(£) = 



1 — I cxp — 



( E - E c f 
2W 2 



(16) 



In order to better characterize the CRF we introduced a new tool, the so-called 
normalized Crab ratio. The Crab ratio is simply the ratio between the source count 
rate spectrum and the count rate spectrum of the Crab Nebula. As this second 
spectrum is known, with great accuracy, to be free of features and to be modeled 



BeppoSAX observation of 4U0115+63: Pulse Phase 0.2-0.3 
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FIGURE 3. Normalized Crab ratio for the BeppoSAX observation of the transient 
AXP 4U0115+63. Both the two high energy instruments, HPGSPC (black marks) 
and PDS (grey marks), are shown. Note the presence of at least four cyclotron 
harmonics |Ga|. 



at first order with a power law in a very broad energy range, this ratio is quite well 
suited to enhance the presence of features in the spectrum. Furthermore the ratio 
is in first approximation independent from the calibration of the instrument. 

In order to enhance the deviations from the continuum we multiply the ratio by 
a E -21 power law, that is the functional form of the Crab Nebula spectrum, and 
we divide by the functional describing the continuum shape of the source (from this 
the name normalized Crab ratio). The procedure is described in Fig. |3J where we 
plot the result of each different step used to obtain the final result in the case of 
4U0115+63. Note the presence of up to four cyclotron harmonics in the spectrum 

By performing the normalized Crab ratio on all the AXPs listed in Table 
it is immediate to observe that higher the CRF energy, broader the feature ||7 
The correlation between E c and the CRF FWHM is quite evident in Fig. [|, and 
is easily understood in terms of Doppler broadening of the electrons responsible of 
the resonance, and holds for all the sources displaying single CRFs (see Eq. g)Q 

2 We assume that CRFs are produced quite close to the NS surface, therefore neglecting the 
gravitational redshift, which shifts the centroid energy by a factor (1 + z) _1 = (1 — 2GM/R.C 2 ). 
See discussion in [pj . 
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FIGURE 4. FWHM vs centroid energy for the cyclotron features observed by Bep- 
poSAX and listed in Table [l]. We include the OSSE measurement on the 1994 out- 
burst of A0535+26 g . Note that the linear correlation has been computed without 
taking into account the A0535+26 point. 



It is important to stress that this relation does not hold in presence of multiple 
harmonics. In other words, it seems that the temperature of the electrons responsible 
of higher CRF harmonics is different from that of the electrons responsible of the 
fundamental CRF. From Fig. and by means of Eq. o we derived that the electron 
temperature responsible for the resonance is in the range ~15-30 keV. This energy 
range is somehow "critical", because in some AXPs (and the effect is particularly 
evident in OAO 1657-40 |lj|) we observe actually two changes of slope in the high 
energy part of their spectra: a first change of slope occurs in the ~ 10-20 keV range, 
while a second steepening occurs for higher energies. This leds to our last issue: 
"anomalous" multiple CRFs. 

There are three sources, namely 4U1907+09, Vela X-l, and GX301-2 that re- 
quire two CRFs in their energy spectra. The anomaly is that (i) the two CRFs 
are not harmonically spaced; (ii) the depth of the "fundamental" is much smaller 
than that of its "harmonics" , (hi) their width does not correlate with their cen- 
troid energy and, more importantly, (iv) there is no trace of the "fundamental" in 
the normalized Crab ratio. Because of this last point we have some doubt about 
the interpretation of them as CRF, expecially because they are all in the critical 
energy range 10-30 keV where we observe the change of slope in the continuum. 
We are therefore inclined to interpret them as due to a not correct modelization of 
the continuum. Another possible explanation could be that they are due to vacuum 



10 



polarization effects (see Fig. fil), but this interpretation requires a more quantita- 
tive analysis. For the three sources discussed above, we did not plot in Fig. || the 
"anomalous" CRF but the one obtained from the normalized Crab ratio. 

4. CONCLUSIONS 

The broad-band capabilities of BeppoSAX have shown that the simple phcnomcno- 
logical spectral laws used to describe AXP spectra in narrow energy ranges are 
inadequate to fit broad-band spectra. The study of AXPs as a class has shown that 
there is a critical region between ^10 and ~30 keV in which we observe a change of 
slope in the continuum. If not well modeled this could give rise to extraneous fea- 
tures that could be interpreted as CRF. Probably a detailed treatment of Compton 
scattering taking into account the effects of the magnetic field could help to solve 
this issue. Also vacuum polarization effects could alter the emergent energy spectra 
and explain the observed "anomalous" CRF harmonics. 

Doppler broadening of the electrons responsible for the CRF is able to explain 
the observed correlation between CRF FWHM and centroid energy, showing that 
the electron temperature is in the range 15-30 keV for all the observed AXPs. This 
correlation does not hold for sources showing multiple CRFs, implying that the 
temperature of the electrons giving rise to higher harmonics could be different. 
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